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We derive analytic expressions of the recursive solutions to the Schrodinger's equation by means 
of a cutoff potential technique for one-dimensional piecewise constant potentials. These solutions 
provide a method for accurately determining the transmission probabilities as well as the wave 
function in both classically accessible region and inaccessible region for any barrier potentials. It is 
also shown that the energy eigenvalues and the wave functions of bound states can be obtained for 
potential-well structures by exploiting this method. Calculational results of illustrative examples 
are shown in order to verify this method for treating barrier and potential-well problems. 
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I. INTRODUCTION 

Quantum mechanical tunneling and quantum well 
problems have attracted more interest, due to recent ad- 
vances in the fabrication of semiconductor layers and the 
development of high-speed and novel devices [l| . In order 
to better understand the physical properties of a device, 
it is important that one can accurately solve and analyze 
the one-dimensional potential problems. A number of 
methods for solving the Schrodinger's equation and for 
locating the bound states and resonances generated by 
one-dimensional potentials have been developed over the 
past decades 0, S S H, Hi ■ Most of them are based on 
the so-called transfer-matrix approach. There are also 
some methods that are cumbersome to implement, for 
instance, the Monte Carlo method 0] and the finite ele- 
ment method (FEM) [1]. 

In contrast to those numerical methods, a limited num- 
ber of exact analytic solutions are available only for sim- 
ple potential structures; analytic solutions are preferred 
due to its simple form and clear interpretation of the 
physics underlying the process. WKB method, as an ap- 
proximate analytic approach, has been widely used but 
it is restricted to slowly varying potential profiles that 
are continuous. Improved methods such as the modified 
conventional WKB (MWKB)[i] and the modified Airy 
functions (MAF) l3| still fail to provide perfect results. 

Tikochinskyfn'l has derived two nonlinear first-order 
equations replacing the second-order linear Schrodinger's 
equation in a cutoff potential method for one-dimensional 
scattering amplitudes. In this paper we solve these non- 
linear first-order equations analytically and obtain recur- 
sive solutions of the Schrodinger's equation for multi-step 
barrier potentials. We intend to find a method provid- 
ing a general analysis of the one-dimensional problem for 
both classically accessible region and inaccessible region. 
In fact, arbitrarily accurate solutions including the scat- 



tering amplitudes and the wave functions are obtainable 
with these recursive solutions for any potential profile by 
dividing it into many segments, since any continuous po- 
tential problem can be recovered as the segments become 
finer and finer. 



Using the recursive solutions, we will also determine 
the wave functions of quasi-bound and bound states. We 
can solve a potential-well problem by evenly uplifting 
the potential function restricted on the potential-well re- 
gion and its surroundings to construct a resonant barrier- 
potential that can readily be handled as a model poten- 
tial. The resonant barrier-potential constitutes a quasi- 
bound-state problem, therefore each sharp tunneling res- 
onance is locatable so that a quasi-bound-state eigen- 
value can be determined which immediately leads to a 
bound-state eigenvalue being sought for the original po- 
tential well. Then, we can determine the eigenfunction 
belonging to the quasi-bound-state eigenvalue by means 
of the recursive solutions. 



In Section II, we review the basic formalism for 
the one-dimensional problems that allows two nonlinear 
first-order equations to replace the linear second-order 
Schrodinger's equation. In section III, we derive ana- 
lytic expressions with recursive coefficients by solving the 
nonlinear equations for the transmission amplitudes, re- 
fiection amplitudes, and the wave functions for piecewise 
constant potentials. In section IV, a quasi-bound poten- 
tial profile which characterizes a resonant tunneling is 
considered as the model potential of a potential well in 
order to treat bound-state problems. In section V and 
VI, we demonstrate the validity of the method by taking 
some examples, and discuss the calculational results. 
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II. BASIC FORMALISM 

The one-dimensional Schrodinger equation for a parti- 
cle with mass m incident upon a potential energy U{z) 
that asymptotically vanishes can be written as 



Tpix) = V{x)lp{x) 



(1) 



where x = zy^2me/h?, fc^ = K/e, and V{x) = U{z)/e 
with a certain unit energy e chosen for convenience. Here, 
z is the position of the particle in the one-dimensional co- 
ordinate, K is the energy of the particle, and h = h/2TT, 
h being Planck's constant. Thus, the coordinate x, the 
wave number k and the potential V{x) are all dimension- 
less. The dimensionless energy E is defined by 



E 



(2) 



Now. we introduce a cutoff potential that was used in 
Refs. [H [H [H, [ii. The cutoff potential is defined by 



V,iy,x) = V{x)e{x-y), 



(3) 



where 9(x — y) is the step function. If V{x) is replaced 
by the cutoff potential, Eq. ([T|) has the formal solution: 



Ak\: 



'''Wciy,x')ipEiy,x')dx'. 

(4) 

Then, the cutoff reflection and cutoff transmission am- 
plitudes for this problem can, respectively, be written as 



and 

TE{y) 



1 



1 



2ikA 



2ikA 



e'''^'v,{y,x')^E{y,x')dx' (5) 



^—ikx 



V,{y,x')^E{y,x')dx' . (6) 



In Ref. nil it has been shown that the problem can be 
stated as two nonlinear first order equations given by 



dREjy) 

dy 



REiy)]' 



(7) 



and 

dTEjy) 
dy 



-iky 



2ik 



-V{y) [' 



RE{y)\TE{y), (8) 



with the boundary conditions, Re{oo) = and Te{oo) = 
1. The Schrodinger's equation, which is a second-order 
differential equation, has been disassembled through the 
cutoff-potential manipulation into two first-order equa- 
tions, Eqs. dH) and ([8|). Integration of Eq. ([8]) immedi- 
ately leads to 



Te{x) = exp 



2ik 



V{y) [e^^y 



''"yREiy)] dy 



(9) 



The reflection amplitude Re and the transmission am- 
plitude Te for the potential V{x) are Re = Re{—oo) 
and Te = Te{—oo) in terms of which the reflection and 
transmission probabilities are, respectively, given as 



Pr = \Re\ 
Pt=^-^\Te\\ 



(10) 
(11) 



where % = \/ E — V^(oo). 



III. THE PIECEWISE-CONSTANT POTENTIAL 
PROBLEM 

We want to find analytic expressions of the transmis- 
sion amplitude, reflection amplitude and the wave func- 
tion for a multi-step potential profile which contains n 
layers of constant potential in the zero-potential environ- 
ment, as shown in Fig.[TJ Let us solve Eq. ([Tj) for Re{x) 
in jth step region (&j-i < x < bj) for V{x) = Vj ^ 0, 



where Vj is constant for j — I, 



Qe{x) = RE(.x)e 
Eq. ([7]) can be rewritten as 



-2ikx 



, n. By defining 



dQE{x) .Vj 



dx 2k 
to which the solution is 



{QE{x)f 



1 



A, 



-2pjX 



(12) 



(13) 



(14) 



with a constant Aj to be determined by a boundary con- 
dition. Here, Pj is defined by 



Pj 



^Vj - fc2 for Vj > fc2, 



(15) 



iy/k^ ~ Vj for Vj < P. 

Then, we obtain the cutoff reflection amplitude in the 
jth step region: 

"2fc2 



Re{x) 



V, 



1 + i2k 



Vi VI 



1 



2ikx 



On the other hand, the solution of Eq. ^ for Vj 
the jth step region is 

REix) = Cj 



(16) 
in 



(17) 



which can be determined by the continuity condition of 
REix) at the boundary x — bj. 

If Vj+i — 0, Re{x) — Cj+i in the region bj < x < bj+i 
and the continuity of Re{x) at the boundary x — bj leads 
to 



A, 



Vj{Cj+ie-^'^''^ + 1) - 2P - i2kpj' 
VjiCj+ie-^'''''^ + 1) - 2fc2 + i2kpj 



(18) 
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FIG. 1: 

A multi-step potential profile consisting of n steps in the 
zero-potential environment. 



From Eqs. ^ and © 



^E{y, y)^A [e'l^y + e-^'^yREiy)] (24) 
and by differentiating Eq. ([4]) in the region y < x one has 
the relation 



dipE{y,x) 

dy 



^^yv{y) 



2ik 



'ipE{y,y)e 



ikx 



(25) 



2ik 



dy 



But, if Vj+i 7^ 0, the boundary condition gives 



A, = 



_ ( D] - Vj+iPj + ikjVj+i -Vj) \ 2p,b 



where 



D] = VjPj+i 



l + A 



1 - Aj+ie-^Pj+ibj 



(19) 



(20) 



Let us take the index j = n + 1 for the region x > bn- 
Since Vn+i — and Pn+i — ik, we obtain An+i = from 
(UHl). Aj is related to Aj+i through Eq. HH) if Vj+i ^ 0, 
while it is related to the constant reflection amplitude 
Cj+i through Eq. ^ if Vj+i = 0. Therefore, all A/s 
for j = 1, • • • , n can be determined recursively from the 
starting value A„+i — 0. 

The cutoff transmission amplitude Te{x) can be ob- 
tained in the jth step region (&j-i < x < bj) from Eq. ([9]) 
by the help of 



Te{x) = TE{b^j)e-K]i[[pj ~ ik){bj - x)] 
1 - A.er'^P^^^ 



1 - yl,e-2pi2: 



(21) 



where 



Eq. (gl) and Eq. ((251) are both the solution of Eq. (U) 
with the cutoff potential T4(y, a;) but have differ- 
ent amplitudes. Their amplitudes differ by a factor 
—e~'^''yV{y)ipE{y,y)/2ikA, thus in the region y < x 



d^E{y,x) _ e '^yyjy) r 



dy 



2ik 



Rsiy)] ipE{y,x). 

(26) 



Integration of Eq. in the region y < x leads to 



tpEix,x) i})e{x,x) 



tPe{x) 4!e{-oo,x) 

/ px —iky 

exp ( - / ^V{y) [e^'^y + e-'^^i?£;(y)] dy 



(27) 



which can be rewritten, with the aid of Eq. as 



ipEix)^A 



— ik 



Te{x) 



^Re{x)\ 



(28) 



TE{hj) = n exp[(pj-ifc)(6j-5j_i)] 

i=]+i 

1 - Aie-^P''" 

1 - A;e-2p!bi-i 



(22) 



Therefore, the analytic expression of the transmission 
amplitude is 



Te = T£;(-oo) = J^exp [(p; - iA;)(6; - 5/_i)] 
1 - Aie-^P''" 



1 - A;e-2P''"- 



(23) 



One can imagine the insertion of an infinitesimally thin 
virtual layer of zero potential at x, and consider the 
wave function iPe{x) to be the superposition of the for- 
ward traveling wave e*'^^ and the backward traveling wave 
Re{x) e~*'^'^ reflected upon the potential on the right side 
of the virtual layer, multiplied by the factor ATe /Te{x) 
which is fathomable to be the correct amplitude for the 
superposed wave function at x. 



Substituting ([21]) and into Eq. ^Q, we write the 
wave function for the energy E in the jth layer {bj-i < 
X < bj) as 



4 



■(pjEix) = AY[exp[{pi - ik){bi ~ bi-i)] ij 



1=1 



1 - Aie-'^P'''' 



X exp [{p, - ik)ix - bj)] ( ^_^ ^_2p^fc^ 



1 - Aic-^P^ 



Re{x)] 



AY[e^v[{pi-ik){bi~bi^i)] (- 



1=1 



1 - A;e-2Pi&i 



lS-[l-A,exp(-2p,6,)] 



[{~ik + pj) ef^" + A, {ik + pj) e-P^' 



(29) 



Eq. (|29p is an analytic expression for the recursive solu- 
tion of the Schrodinger equation. It is simple and fast 
to calculate Aj^s for any multi-step potential. Also, any 
continuous potential problem may be considered to be 
a multi-step potential problem by dividing the non-zero 
potential region into a lot of small segments such that 
the potential in each segment can be approximated to a 
specific value, for instance, the average value, of the con- 
tinuous potential within the segment. Thus, the analytic 
expression of the wave function given by Eq. (p9)) can 
be used for obtaining an arbitrarily accurate solution for 
any one-dimensional potential. 

In passing, let us check that Eq. ([29]) is consistent 
with the traveling plane waves which are expected in the 
regions x < bo and x > bn, where V{x) — 0. With- 
out loss of generality, one may take the first step region 
{bo < X < bi) such that Vi — 0. Then, the wave function 
in the region x < bi becomes 



i;iE{x) = lim A(e"=" + e-''="i?£(x)) 



X exp[(pi - ik){x - bo)] 



1 - Aie-^Pi 



1 - Aie-2pi6o 
=A [e''^^ + Re e'^'^^^j , (30) 

since Ai — — (Fi/4fc^)i?_E from Eq. ([TBI) for the limit 
Vi 0. On the other hand, if we take the last step 
region 6„_i < x < bn such that V{x) ^ Vn ~ 0, the 
wavefunction in the region x > becomes 



(31) 



since Re{x) = and An — — (T4/4fc^)i?£;(x) from 
Eqs. (HH) and ^ in the Hmit Vn 0. Eqs. ^ and 
([3T|l show that Eq. ([29|l is consistent as the solution of 
Eq. ID) in the zero-potential region. 



IV. QUASI-BOUND STATES 

We consider a resonant barrier structure that con- 
sists of n layers of piecewise constant potential. Reso- 
nances are supposed to occur at energies less than both 



Vi and Vn which are potential values of the end steps so 
that quasi-bound states are involved. The wave function 
ipjE^x) in jth layer for this structure can be obtained 
with the calculated Aj's which are dependent on E. An, 
in the last layer, can be evaluated from with cutoff 
reflection amplitude Cn+i = 0. Inserting the evaluated 
An into (|29p . we find in the last layer for E <Vn 



ipnE{x) oc {-ik + pn)eP"'' 

Vn - 2fc2 - 2ikpn 



(32) 



Vn - 2fc2 + 2ikpn 



{ik + Pn 



where the first term exponentially increases with increas- 
ing X, while the second term exponentially decreases. 
One can see that if p„(6„ — 5„_i) is large enough, the 
first term is negligible in magnitude compared with the 
second term, so the wave function exponentially damps 
with increasing x in the region bn-i < x < bn- 

Now, let us consider the wave function in the first layer 
which is given from (|29p by 



tpiE{x) = A 



2ike''''"> 



Vi{eP^'"> - ^le-Pi^o) 
[{-ik + pi)eP''' + {ik + pi)Aie-P'' 



(33) 



where Ai is obtained through the step potential profile 
of the resonant barrier structure. The wave function 
given by Eq. p3p exponentially increases with increasing 
X when — 0; Ai vanishes at the energies of quasi- 
bound states. 1 = may be used as the quantization 
condition for the energy levels of a particle confined in 
a potential well. In the practical calculations, transmis- 
sion resonances occur at the energies where \Ai\ takes a 
minimum value in a sudden dip, if bi — bo and 6„ — bn~i 
are large enough. Such resonances indicate quasi-bound 
states. 

Now, a potential well can be treated by uplifting the 
potential function restricted on the region of the well and 
its surroundings to form a resonant barrier structure that 
can be handled in our approach. This resonant struc- 
ture is used as a model potential of the potential well 
under consideration. Then, each resonance occurred in 
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FIG. 3: 

Resonances for a model potential of harmonic oscillator 
given by Eq. (|34|l . Each sharp dips on the solid line of the 
curve In ([ylil) exactly agrees with a resonant peak on the 
dashed line of the curve In (Pt). 



FIG. 2: 

(a) Transmission probabilities of double-barrier potentials: 
the solid line and the dashed line are for the square and the 
gaussian double-barrier potential, respectively. Their 
potential profiles are shown in the inset graph, (b) The 
curve of In for the double square-barrier potential. Its 

profile is in the inset graph. Each dip of the ln(|j4i|) curve 
below the barrier exactly agrees with the resonant peak of 
transmission probability shown in (a). 

this model-potential problem corresponds to an energy 
eigenvalue for the potential well, and the real part or 
the imaginary part of the quasi-bound-state function ob- 
tained for a resonant energy is the eigenfunction of the 
bound state belonging to the corresponding eigenvalue 
for the potential well within a normalization constant. 



V. SOME EXAMPLES 

In order to demonstrate the validity of our method, we 
consider quantum barriers and quantum well potentials. 
Here, we examine the derived recursive expressions of 
the transmission amplitude and the state wave-function 
given by Eqs. and (pS)) . respectively, by calculating 
Aj^s for several potentials. 

First, we consider rectangular and gaussian double- 
barrier potentials which are shown in the inset of 
Fig. [2Ia). Here, the unit energy e is chosen to be the 
maximum potential energy of the barrier. Since the rect- 
angular double barrier structure consists of three layers of 
constant potential, the analytic expression of the recur- 
sive solution ijjjEix) for this structure can be obtained 
with calculated Aj^s which are dependent on E. One 
can determine A3 from Eq. ([T5)l with the cutoff reflec- 



tion amplitude C4 = Re{x) = in the region a; > 63, 
and determine C2 for the step region of j = 2 where 
V2 = by evaluating Re{x) at o; = &2 in the step region 
of j = 3 through Eq. (fT6|) . Then, Ai can be also ob- 
tained through Eq. (flS]) using the value of C2 . Fig. [2jb) 
shows the curve of ln(|v4i|) for the rectangular double 
barrier structure. Each dip of the curve coincides with a 
transmission resonance for the energy below the barrier 
{k < 1) as shown in Fig. [S] Such a sharp dip of Ai indi- 
cates that an exponentially increasing term of the wave 
function is dominant in the first barrier region due to a 
neghgibly small \Ai\ as was considered below Eq. 
so that the transmission probability becomes large in the 
narrow energy region. The resultant transmission proba- 
bilities for the two barrier profiles have similar resonance 
patterns as shown in Fig. [2ja). We have confirmed that 
the solid line for the rectangular double barrier perfectly 
agrees with the result obtained by the traditional method 
for solving the Schrodinger's equation for step potential 
problems with boundary conditions imposed. 

Next, we consider a potential model by which one can 
obtain eigenvalues and eigenfunctions for a harmonic os- 
cillator. We take the unit energy to be e = hjjj/2, where 
Lo is the characteristic angular frequency of the harmonic 
oscillator. Then, the dimensionless potential in Eq. ^ 
is Vho — for the harmonic oscillator. In order to apply 
our method to the harmonic oscillator in an energy range 
< -B < 16, we employ a model potential: 

{x^ for |a;| < 4 
16 for 4 < < 10 (34) 
for \x\ > 10, 

as shown in the inset of Fig. [3l This model has a dou- 
ble barrier structure and represents a quasi-bound prob- 
lem in which the transmission of a particle through the 
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FIG. 4: 

The wave function of the lowest quasi-bound state {E = 1.0) 
for the potential given by Eq. (|34p for three regions with 
different units. The solid line and dashed line are the real 
and imaginary parts of the wave function, respectively. 




model potential reaches its peak when the energy is res- 
onant with one of the energy levels of the quasi-bound 
states in the region between the two potential barri- 
ers. The potential in Eq. (|34p is expected to involve low 
quasi-bound-state energies that approximately equal the 
bound state energies of the harmonic oscillator, respec- 
tively. Each sharp dip of ln(|Ai|) suppresses the term 
exponentially decreasing with increasing x in the wave 
function in Eq. (|29p allowing the exponentially increas- 
ing term to be dominant in the first barrier, and thus 
occurs at the same energy as a resonance peak which lo- 
cates a quasi-bound state, as shown in Fig. [3l Energies 
of the quasi-bound states have been determined very ac- 
curately to fit in the exact eigenvalues of the harmonic 
oscillator £"„ = 2n + 1 with rt = 0, 1, 2, • • • for the low 
lying states. However, there are very small errors for the 
upper levels, for instance, n = 5,6, 7, because the model 
potential in Eq. is different from the potential of the 
harmonic oscillator for the region > 4. 

Once a quasi-bound-state eigenvalue is determined, the 
wave function of the state can be obtained from Eq. (^5)1 . 
The lowest quasi-bound state for the model potential in 
Eq. (pi]) corresponds to the ground state for the har- 
monic oscillator. The wave function of the state is sepa- 
rately shown in Fig. [4] for the regions, —30 < a; < —10, 
— 10 < a; < 10, and 10 < a: < 30 with different units. 
The solid line and dashed line are, respectively, the real 
part and the imaginary part of the wave function. The 
forward and reflected plane waves are superposed in the 
region x < —10, and only the transmitted forward plane- 
wave propagates in the region x > 10. One can see in 
Fig- mb) that the wave function exponentially increases 



FIG. 5: 

The wave functions of three quasi-bound states of energies, 
(a) E = 3.0, (b) E = 5.0 and (c) E = 7.0, for the potential 
given by Eq. (|34|l . The solid lines and dashed lines are the 
real and imaginary parts of the wave functions, respectively. 



with increasing x in the first barrier region at the reso- 
nance energy but exponentially decreases in the second 
barrier region which is sufficiently wide as explained ear- 
lier, forming the lowest quasi-bound-state wave function 
in the binding potential region between the two barriers. 
Wave functions of the next three quasi-bound states for 
the same model potential that correspond to the three 
lowest excited states in the harmonic oscillator problem 
are also shown in Fig.O Either the real part or the imag- 
inary part of a quasi-bound-state wave function can be 
taken as the wave function of the corresponding bound 
state for the harmonic oscillator. 

For our final example, we consider a symmetric double- 
well potential given by 



Vdw{x) 



x'^{x^ 




16) for |a:| < 4, 
for \x\ > 4. 



(35) 



For treating this problem with our methods, we take a 
model potential of Vdw by uplifting the potential function 
on the region < 6 by 64: 



x'^{x^ 
Vdw{x) = { 64 




16) + 64 for |a;| < 4, 

for 4 < |a;| < 6, (36) 
for \x\ > 6, 
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FIG. 6: 

Logarithmic transmission probabilities for double-well model 
potentials given by (a) Eq. (|36|l and (b) Eq. (|37p . Each peak 
represents transmission resonance. 

which is shown in the inset of Fig. [6{a). The local max- 
imum of the model potential located at a; = is 64, 
and the two local minima are located at x — 2\/2- The 
calculated curve for In Pt is plotted in Fig. [IJa). As is 
well-known, each peak comprises a pair of states, an even 
state and an odd state, whose energies are the very same 
especially for a peak of lower energy. But they are not 
perfectly degenerate, since a particle in one potential well 
can tunnel into the other well through the middle bar- 
rier which is finite. The energy difference of two states 
comprised in a peak increases with increasing energy, and 
thus the seventh pair of states which are just below the 
barrier appears to be separated. 

In order to confirm that two states are comprised in a 
peak, we add a small asymmetric potential function to 
the potential in Eq. 



Vdw{x) + V'{x) = Vdw{x) +tSinhx+l. (37) 



Using this potential, we obtain a curve of In \Pt\ in which 
there are seven pairs of quasi-bound states as shown in 
Fig. 6(b). Each pair is split into two peaks due to the per- 
turbing potential V' (x) from a peak containing two states 
for the double-well potential Vdw shown in Fig. 6(a). 
The lower state and higher state in each pair are expected 
to be confined in the left well and right well, respectively. 
In order to find bound-state wave functions for the dou- 



FIG. 7: 

The wave functions of three bound states for the double well 
potential, V{x) = x^{x^ — 16) -|-tanhx; (a) the ground state, 

(b) the eighth excited state, and (c) the ninth excited state. 
These wave functions have been obtained by using the 

model potential given by Eq. (|37|l as explained in the text. 



ble well potential given as V{x) = x^(a;^ — 16) -I- tanhx, 
one may take either the real parts or imaginary parts 
of quasi-bound-state wave functions determined from 
Eq. (p9)) for an appropriate model potential such as 
Eq. p7p . The real parts (or imaginary parts) can be con- 
sidered as the bound-state eigenfunctions within a nor- 
malization constant. Real parts of the wave function of 
the lowest {n — 0), the ninth (n = 8), and the tenth 
{n ~ 9) quasi-bound states, which are taken to be the 
wave functions of corresponding bound states, are plot- 
ted in Fig. 7(a), (b) and (c), respectively. Fig. 7 shows 
that the particle which occupies the lowest state or the 
ninth state stays in the left potential well, while the par- 
ticle occupying the tenth state stays in the right potential 
well, as expected. 



VI. DISCUSSIONS AND CONCLUSION 

Analytically solving two nonlinear first-order differen- 
tial equations equivalent to the Schrodinger's equation, 
we have obtained the recursive solutions of the cutoff 
reflection amplitude Re{x), the cutoff transmission am- 
plitude Te{x) and the state wave function ipEi^) for a 
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multi-step potential. One can use the recursive solu- 
tions in Eqs. (p3| and ((29)l for any piecewise constant 
potential to obtain transmission amplitudes and eigen- 
functions, respectively. If an approximate multi-step po- 
tential consisting of sufficiently many layers is substituted 
for a smooth potential profile, calculations with the re- 
cursive solutions are accurate and fast. 

For treating a potential well, we take a model potential 
with two side barriers that contains a partial profile iden- 
tical to the profile of the potential well in the region be- 
tween the side barriers. In the model potential problem, 
resonance occurs at the energy of a quasi-bound state, 
so the energies can be determined by locating the trans- 
mission resonances. Once a quasi-bound-state energy is 
known, the eigenfunction for the energy can also be de- 
termined. However, the real part and the imaginary part 
of the wave function of a quasi-bound state exponentially 
increases at a same rate in the first barrier and exponen- 
tially decrease in the last barrier as explained earlier, and 
are the same within a factor tan <j) in the region between 
two side barriers, 4> being a phase angle. Therefore, ei- 
ther the real or the imaginary part can be taken for the 
wave function of the corresponding bound state with a 
normalization constant. 

We have considered several examples to demonstrate 
the validity of our method with the recursive solutions. 
First, we have calculated the transmission probabilities 
of a particle incident upon a rectangular and a gaussian 
double-barrier potential which show typical features of 
resonant tunneling in Fig. [D^a) . The rectangular double 
barrier simply consists of three layers, the left barrier be- 
ing the first potential step, and the feature of resonant 
tunneling can be explained by the calculational results 
of Ai in which information of the barrier structure is 
included, li Ai = for an energy lower than the barri- 
ers, the wave function does not contain exponentially de- 
creasing term with increasing x in the first barrier region 
so that it exponentially increases in the region, and thus 
the transmission probability reaches maximum at the en- 
ergy. Practically, the sharp dips of the ln|Ai| curve in 
Fig. 2(b) represent the resonances associated with quasi- 
bound states in the rectangular well between the barriers. 

The harmonic oscillator problem can be dealt with 
by solving a quasi-bound-state problem with the model 
potential given by Eq. ([M)) . The calculational results 
for In 1 and In Pt are shown in Fig. [Sj so the quasi- 
bound-state energies which may be regarded as the en- 
ergy eigenvalues of the harmonic oscillator can readily 
be determined by locating the resonance peaks of In Pt 
or the sharp dips of ln|Ai|. The resultant peaks per- 
fectly agree with the exact eigenvalues for low energy 
states, but deviate a little from the exact values for high 
energy states; Eq = 1.000, Ei = 3.000, E2 = 5.000, 
E3 = 6.9999 E4, = 8.999, E5 = 10.994, Eq = 12.970, and 
El = 14.857. Since the deviation originates from the 
approximated model potential which does not agree with 



the harmonic-oscillator potential on the high potential re- 
gion I a; I > 4, one may get more accurate eigenvalues for 
higher states by using a larger model potential which co- 
incides with the harmonic-oscillator potential in a wider 
region. 

We have applied our method to a symmetric double- 
well potential given by Eq. ([35|) . using Eq. ([36]) as its 
model potential. It is well-known that an even state and 
an odd state are folded together belonging closely to an 
eigenvalue if they are deeply bound in a symmetric dou- 
ble well whose middle barrier is strong enough so that the 
wave functions nearly vanish in the barrier region. The 
even state and odd state belonging to the same eigen- 
value can be combined to constitute two bound states 
that are localized in the left well and right well with the 
same eigenvalue, respectively. Thus, one may consider 
each resonance peak to correspond to a pair of states 
confined in one and the other well, as well as a pair of 
even and odd states. When an asymmetric perturbing 
potential is added to the double-well potential, each res- 
onance peak has been split into two peaks, as shown in 
Fig.[6fb). Since the potential is not symmetric due to the 
perturbation, even and odd states can not be sustained. 
Therefore, each resonance peak in Fig. represents 
a state confined in one of the two wells. As expected. 
Fig. [7] shows that among the two states separated from 
each other due to the perturbation, the lower energy state 
is confined in the left well which is a little deeper than 
the right well, while the higher energy state is confined in 
the right-well region. The ninth (n=8) and tenth (n=9) 
states are the fifth states in the left well and the right 
well, respectively, so the numbers of nodes in their wave 
functions are identically four as in Fig 7 (b) and (c). 

Here, more accurate determination of quasi-bound- 
state energies is required to obtain eigenfunctions of the 
lower quasi-bound-state confined in the right well, be- 
cause the resonance is sharper for the particle tunneling 
into a deeper quasi-bound state or a state confined in the 
right well through the stronger barrier or two barriers 
(the first barrier and the middle barrier); for instance, 
the determined eigenvalues Eq = 5.601849104, Ei = 
7.58342367856952, Es = 46.290706 and Eg = 48.1538536 
should be used in obtaining the wave functions, where Ei 
is required to be the most accurate. Here, each bound- 
state energy for V{x) = x^(a;^ — 16) -I- tanhx is the cor- 
responding quasi-bound-state energy for V[)^ minus 65. 
Unlike the model potential of the harmonic oscillator, 
no approximation is involved in the calculation with the 
model potential of this double well potential. 

Although here we have presented results for double 
barriers, a harmonic oscillator, and double- well poten- 
tials, the analysis can easily be applied for accurate cal- 
culations to general potential barrier and quantum well 
structures. The analysis could be extended to the one- 
dimensional problem of a particle that has a coordinate- 
dependent mass, by imposing the boundary conditions 
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on the step potential structure. In addition, the three- 
dimensional Schrodinger equation for a spherically sym- 
metric potential can also be solved using the present 
method because such a three-dimensional problem can 
be reduced to a one-dimensional problem. 

In conclusion, we have described a method for solv- 
ing the one-dimensional Schrodinger equation by means 
of the analytic expressions of recursive solutions derived 
for piecewise-constant potentials. The recursive solutions 
provide a general analysis of one-dimensional scattering, 
quasi-bound and bound state problems. We have demon- 
strated the validity of the method by taking some exam- 
ples to show consistent and predictable calculational re- 
sults. The calculations with the recursive solutions were 
rapid and accurate even for smoothly varying potentials. 
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